Upper critical field in a trigonal unconventional superconductor UPt3 
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A theory of the upper critical field in a trigonal superconductor with a two-component order- 
parameter is developed. We demonstrate the existence of sixfold modulations of the upper critical 
field in the basal plane. The form of the angular dependence of H 2 in the plane of the main 
symmetry axis is studied in the whole range of phenomenological rigidity coefficients. A qualitative 
dependence of these coefficients on the form of the Fermi surface is established. The results obtained 
are discussed in connection with the measurements of the direction dependence of H C 2 in the heavy 
fermion superconductor UPt3 that has been recently found to possess trigonal crystalline structure. 

PACS numbers: 74.20.De, 74.70.Tx 



I. INTRODUCTION 

The heavy-fermion superconductor UPt3 is tradition- 
ally held to be hexagonal. Based on high-energy x- 
ray diffraction and transmission electron, microscopy, the 
crystal structure of UPt3 have recently^ been argued to 
be determined to be trigonal. 

Trigonal deviations from the "ideal" hexagonal symme- 
try are of the order of 1%. That is why one can expect 
only small changes in the quantitative interpretation^ of 
the numerous measurements of the thermodynamic and 
kinetic values based on the Em two-dimensional super- 
conducting state. The same one can say in respectppf the 
interpretation-of the flux-line lattice experimental pro- 
posed in Ref.cl. It seems to be interesting however to 
concentrate efforts on the study of the properties that 
can demonstrate clear indubitable distinction between 
the two component superconductivity in hexagonal and 
trigonal crystals. The simplest example of such a prop- 
erty is orientational dependence of the upper critical field. 

In the present article we develop Ginzburg-Landau the- 
ory of the orientational dependence of the upper criti- 
cal field for a two-component superconducting state in 
a crystal with a trigonal symmetry. We demonstrate, 
in particular, that the upper critical field in a trigonal 
unconventional superconductor exhibits six-fold modula- 
tions in the basal plane, this conclusion remaining intact 
in the presence of an orthorhombic anisotropy in the di- 
rection fixed parallel or perpendicular to the magnetic 
field. 

It is interesting to note that such small modulations 
(in experiments \5H^ 2 ex / H C 2\ < 0.02) have been ob- 
served in UPt« by angular resolved magnetoresistance 
measurements! 

This experimental observation then encountered diffi- 
culty in explaining in the framework of any model of the 
superconducting phases in a hexagonal superconductor. 
BurlachkovO showed that is isotropic for fields in the 
basal plane near T c if the order parameter rj transforms 
according to the two-dimensional irreducible representa- 



tion Ei of the hexagonal point group D§h- He pointed 
out that is also isotropic for E 2 . 

The hexagonal anisotropy in the GL functional then 
shows up only in the terms of the sixth order in gradi- 
ents, giving that the resulting hexagonal anisotropy of 
H^2 vanishes as (1 — T/T c ) 3 near T c just as it does in a 
hexagonal superconductor described by a one-component 
order parameter. _ _ 

Shortly after the experimentsQ, Minee\fl demonstrated 
that in the model of two accidentally nearly degenerate 
ID representations the hexagonal anisotropy of also 
vanishes, as (1 - T/T c ) 3 near T c . 

Saulsil considered the hexagonal oscillations in the 
model of a two-dimensional representation for the super- 
conducting order parameter coupled to an in-plane an- 
tiferromagnetic (AFM) order parameter m s . Large Zee- 
man energy oriented m s near perpendicular to the ex- 
ternal magnetic field. Hexagonal anisotropy of H c2 then 
arose from the weak in-plane anisotropy energy of the 
AFM state through the coupling of the superconducting 
order parameter with m s . 

On the other hand, the existence of the six-fold 
anisotropy oc (1 — T/T c ) of the upper critical field in 
the basal plane of a trigonal crystal is quite natural, as 
we will see, and does not require any additional assump- 
tions. It has been noticed for the first time numerically 
in the theoretical papero, the main goal of which was 
to study another type of anisotropy of the upper critical 
field — in the plane passing through the high-symmetry 
axis of an exotic trigonal superconducting monocrystal 
of Cui. 8 Mo 6 S 8 . 

We also find the change of the form of H C 2 in the plane 
of the main symmetry axis with the change of the pa- 
rameters of the GL functional placing emphasis on the 
role of trigonality. We determine as well the qualitative 
dependence of the values of these parameters in the mi- 
croscopic weak-coupling BCS theory on hexagonal and 
trigonal contributions to the shape of the Fermi surface. 

The paper is organized as follows. In Sec. II we write 
out phenomenological GL differential equations on the 
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spatial form of the order parameter used to find upper 
critical field in the basal plane in Sec. fill and in the 



plane of the main symmetry axis in Sec. [V. The rigidity 
coefficients in the gradient energy of the GL functional 
are calculated microscopically in the weak-coupling BCS 
theory as functions of the initial trigonal anisotropics of 
the Fermi surface and the order parameter in Sec. [v|. 
Relevance of the theory to the experimental data known 
to date is discussed in the Discussions section [yj. 



II. GL EQUATIONS 

In a trigonal crystal the GL functional up to terms of 
the second-order in the gradients and in the order pa- 
rameter 77 = (771,772) has the form 

T = a\rj\ 2 + kiIAi&I 2 + K 2 \D i T li \ 2 + K 3 (A%)(£>j?7i)* 

+ K 4 \D z f] t \ 2 + K 5 ((A^,77 3 )(L> z 77fc)* +C.C.) , (1) 

where 



a = a (T - T c )/T c , 
the covariant space derivatives 

D = -id + (2tt/$o)A 



(2) 



(3) 



substitute the ordinary ones in the presence of the mag- 
netic field H = d x A, a k is a vector (<r z , —a x ) of the 
two Pauli matrices. <&o/2tt = hc/2e is the flux quantum. 

The condition of the positive definiteness of the func- 
tional (|j) requires that—the following constraints on the 
coefficients are met (cfji3) 

«1 > |« 3 j, K123 > 0, K 4 > 0, Kl23 + «2>0, 

2k\ < K13K4, (4) 

where, e.g., K13 denotes K\ + K3, 

The GL equation ST/Sr]* — is obtained by variation 
of the energy functional with respect to the order param- 
eter. Before varying the functional ([j]) we note that the 
vector (<J Z ,— <j x ) of the two Pauli matrices has the prop- 
erty ofj?7fc = a\\J}k which can be checked directly. Using 
this the term proportional to K5 in the GL functional can 
be rewritten as 

«6 ((A4%)(£>**>)* + (D^a^DirjkT) . (5) 
Performing variation now we get 

-at] = (K 1 D 2 + n i D 2 )r 1 +{K2DD l + K 3 ,D l D)ri l 
+ K5 (D z Di + DiD z ) a-ijrij. (6) 

III. H C 2 IN THE BASAL PLANE 

For the magnetic field directed along the unit vector 
n = (cos ip, sin ip, 0) in the basal plane it is convenient 



to choose the gauge of the vector potential in the form 
A = —Hzt, where f = z X n. Then we may seek for 
solutions of equation (|^) dependent only on z and hence 
put D = (27t/$ )A and D z = -id z in (|). Measuring 
distances in units of the length y / $o/27rff and denoting 

A = -a$ /27riJ, (7) 

we arrive at: 

AT7 = (-Kid 2 + KlZ 2 ) T) + K 23 Z 2 T (TTj) 

+ iK 5 -ficr i j(l + 2zd z )'f]j. (8) 

Rotating vector 77 in the basal plane by the angle ip, i.e. 
changing variables from 771,772 to 7711 = hr),rj± = TT7, we 
finally get the system of two second order linear ordinary 
differential equations 



AT7 = —Kid z r\ 



K! 

K123 



o 

Z 1} 



IK 5 



SmSip -C0s3</?\ , „\ /n\ 

o ■ o (1 + 2zd z )T]. (9) 
cos3^ sm3(y9 I ' ' w 



In a hexagonal superconductor K5 = and equations 
on 77|| and ?7_i_ decouple, each being a Schrodinger equation 
of a harmonic oscillator with its own rigidity coefficient 



( — Kid 2 + K\Z 2 } 

-Kid 2 + K123Z 2 



V-L 



= A||77||, 
= A_L77 ± . 



(10) 

(11) 



Its eigenvalues with the boundary condition of vanishing 
at infinity are 



A|| = y / Ki7t4(2n + 1), 
A_l = ^Ki 23 Ki(2n + 1). 



(12) 
(13) 



When lowering H at a fixed temperature T < T c (which 
corresponds to raising A), the solution realizes for the 
first time at the lowest lying level X n at n — 0. I.e., 
a nucleus of the superconducting phase appears for the 
first time at 



H C 2\\ — 



$ a Q T c -T 



2n y/KlKi T c 



(14) 



if K23 > (and then it is a phase with the order parameter 
t] parallel to the direction of the applied field) , or at 



H,, 



$0 



o: 



T r —T 



2n y/Ki23Ki T c 



(15) 



if K23 < 0, and then it is a nucleus of the phase with 
t] _L H. In any case, neither ( |l4| ) nor (|l5| ) depends on 
the direction of the field, so the upper, critical field is 
rotationally invariant in the basal planeO. 

In a trigonal crystal the coefficient K5 is non-zero, and 
the two equations in (^) cannot be decoupled by a rota- 
tion of the coordinate axes. The influence of the term 
proportional to k 5 can be taken into account as pertur- 
bation. 



Unfortunately, this does nor lead to a concise analytical 
expression. The reason is that in first order of perturba- 
tion theory the correction to the ground state eigenvalue 
of the system (||) is given by the averaging of the pertur- 
bation o perato r oc K5 over the u nperturb ed ground state 
(ipo(z, ^/ki/k 4 ), 0) or (0, ip (z, ^m/^)) depending on 
the sign of K23. Here 



(16) 



is the ground state of the Schrodinger Eq. (ttOh or (|11[) . 



The average of the operator l + 2zd z over the state ( |16| ) is 
zero. So one need to account for second order corrections, 
which involve an infinite sum of matrix elements of the 
operator 1 + 2zd z between zeroth ([TJ]) and higher lying 
states. 

One can, however, accomplish perturbational calcula- 
tions with respect to the parameter K23 , with K5 generally 
staying arbitrary. Then one can obtain expressions for 
the limit K5, K23 -C ki (see Table I for values of the coef- 
ficients K23, k 5 in a BCS theory). 

One proceeds by noting that in either of the cases 
K23 = or k 5 = by some rotation of the coordinate 
axes the two equations (||) decouple. So in both cases 
the upper critical field corresponds to the lower of the 
eigenvalues of the ground states of the two equations. 

Rotating the coordinate system of Eq. (||) by — 7r/4 — 
3(^/2, one arrives at 



1 

-1 



«23 2 

2 



1 + sin 3^5 — cos 3ip 
— cos 3ip 1 — sin 3ip 



1 + 2zd z )r] 
V- (17) 



In the zeroth order in K23 the two equations decouple 
and with the help of the substitution 

m (z) = e** 2 771(2), 77 2 (z) = e 2 m ( z ) (i 8 ) 
reduce both to 

(-K 4 dl + k 1 R 2 z 2 ) rj = A77, (19) 
where we introduced the notation 



R = yjl - k\Ik x k 4 . (20) 
So that the ground state with the eigenvalue 



t/kikIR = J K1K4 — Kg (21) 
is twice degenerate with the wave functions 



e l «4 Z 2 i )0 (z,R x ff),0) , (0,e Ca fo( z , R /£) 
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azimuthal angle cp, rad 

Figure 1: Upper critical field in the basal plane of a two- 
component trigonal superconductor for several values of the 
coefficients K3 and K5. See explanations in the text. 



Using the first order theory of perturbation of a degen- 
erate level, we find 



H, 



R 



K23 1 ± \J sin 2 Zip + R 3 cos 2 3tp 



4ki 



R? 



(23) 



where one should take sign plus for K23 < and sign 
minus for K23 > 0. 

We recall that this formula is valid for small K23 <C ki 
and arbitrary K5. For small K5 and small K23 we get by 
expansion that the first term dependent on ip is 



\ K 23\ 



■ cos 3p 



(24) 



(22) 



while the constant part of H C 2 acquires additional terms 
0(k23/ki) and 0(n 2 / K4K4). 

Eq. (|2J) shows that six-fold modulations of the upper 
critical field in a trigonal crystal have very small ampli- 
tude for small K23 and K5, while their maxima occur at 
(p = (2n+l)7r/6 (i.e., at directions of H || a*) irrespective 
of the signs of k 2 3 and k 5 . 

Numerical calculations confirm these conclusions. 
They also show that the dependence H&fa) rests 



smooth, "sinusoidal" for almost all values of the rigid- 
ity coefficients except for both n 2 3 and K5 large — each 
close to the upper limit of its domain of definition (Q), 
when the maxima of H&fa) at ip — (2n + l)ir/6 become 
sharp. 

On Fig. [l] we presented several curves H C 2(f) obtained 
numerically for K4 = 1.67ki and several K23, K5: 3 groups 
of solid lines show curves for k 5 = 0.3ki, 0.6«i, and 0.9ki 
respectively from bottom to top. Inside each group 4 
curves correspond to K3 = 0, 0.2ki, 0.5ki, and 0.999ki 
(from the upper line down to the lower). Two dashed 
lines display H C 2(tp) for K3 = — 0.15«i and K5 = 0.3ki 
(the lower of the dashed curves) and 0.6ki (the upper). 

An alternative variational method of finding sinusoidal 
corrections ( |2~^ ) to the upper critical field, which also pro- 
vides some insight into the behavior of H C 2 in a variant 
of the GL theory with an orthorhombic symmetry term, 
is analyzed in Appendix A. 



IV. H c2 IN THE PLANE OF THE MAIN 
SYMMETRY AXIS 

Upper critical field in the ac plane rfor a-trigonal su- 
perconductor was considered in paperscl andlla in connec- 
tion with the Chevrel phase compound Cui.gMoeSs- We 
would like to emphasize here the role of the coefficient 

In the case of a magnetic field lying in the etc plane the 
GL differential equation on r\ analogous to (||) is more 
complicated. It can be somewhat simplified if one as- 
sumes 



K 2 = K 3 . 



(25) 



This equality comes out in the microscopic calculation of 
the rigidity coefficients (see Sec. |V ), it was also presup- 
posed in bothtl andEI Conditions (41) become then 



Kl > 0, K4 > 0, k 5 < 



< k 3 < ki. (26) 



However, even with k 2 = K3 the GL equations allow 
analytical solution only for k 2 — K3 — 0, K5 = 0, where 
the upper critical field 



H. 



(0) 



a 



2tt 



ki (ki cos 2 9 + K4 sin 2 



(27) 



turns out to be the same as in the case of a one- 
component superconductivity in uniaxial crystals de- 
scribed by the effective mass model (see, e.gEI). 

The ground state of the Schrodinger equation whose 
eigenvalue is given by ( p7|) is twice degenerate. And the 
influence of K3 and K5 can be accounted for by perturba- 
tion theory of degpaerate levels. Corresponding theory 
was constructed ir£3. The formulas, though get more ac- 
curate in the second order taking into account crossing 
of the first-order levels, are barely citable. 




jt/2 

polar angle 0, rad 



Figure 2: Upper critical field H C 2 in a hexagonal superconduc- 
tor with two-dimensional order parameter for a magnetic field 
lying in the plane of the main symmetry axis as a function of 
the angle 9 between this axis and the magnetic field. Curves 
in the two sets correspond to equally spaced values of K3 run- 
ning the interval from to 0.999ki (upper plot, from top to 
bottom) and from to -0.333fi;i (lower plot, from bottom to 
top). 



Unexpectedly, a simple first-order expression 



AO) 



1 + — + -rS sin 29 

K\ 2k 



l + ^ CO s 2 0-^sin20 
K 2k 



(28) 



where is given by (p?j) and 



K = Ki cos 2 9 + K4 sin 2 9, 



(29) 



is in no much worse accord with the values obtained nu- 
merically than the bulkiest of those ir£3. 

Note that variational fit with the trial function r/ = 
(Vi: 7 l2)e~ Az / 2 with three free parameters givestl an ex- 
pression which in the domain of applicability of pertur- 
bation theory coincides with (pg|). 

The dependence of the upper critical field in a hexag- 
onal crystal (i.e., for K5 = 0) for a magnetic field lying 
in the ac plane on the angle 9 between the direction of 
magnetic field and the c-axis is depicted on Fig. || for 
several values of the coefficient K3 from its definition in- 
terval (26). K4/K1 was put equal to 1.67. Only the range 
of angles from to ir is shown because reversing direction 
of the magnetic field, which corresponds to the substitu- 
tion 9 — » 9 + ir, does not change the value of H c2 resulting 
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7l/2 71 

polar angle 6, rad 



Figure 3: H C 2 in a trigonal crystal («g 7^ 0) for K3 = 0.4ki 
(upper plot), K3 = (middle plot) and kz = — 0.15fi;i (lower 
plot). In each set of the curves the coefficient K5 runs over an 
equally spaced set of values in the interval from (lower curve 
in the set) to the upper boundary y' K13K4/2 of the possible 
values of K5 (upper curve in the plot). 

from the functional ([j]) with (0) , meaning that the curve 
on Fig. § should be continued periodically. 

The curves in each set correspond to the equally spaced 
values of the phenomenological coefficient K3 running the 
range between K3 = and K3 = 0.999ki (upper plot), 
and between K3 = and K3 = — 0.333ki (lower plot) 
inclusively. 

The behavior of H c2 with the increase of is de- 
scribed (for large | K3 1 of course only qualitatively) by 
( p8|) with K5 put equal to zero. Since 

1 / COS 2 8 „ Ha „ . , 

— / — — = 1 + — tan 2 (9 > 1, 30 

K\ I K K\ 

we conclude that for K3 < 0, when the first term in ( |2g|) 
is the larger, the curve of the dependence of if C 2 on 8 is 
scaled with the increase of \k^\ roughly as 

1/(1 + k 3 /ki) (31) 

without changing its overall shape. Note also that 
H c2 {0 = f ) then equals @. 




Figure 4: Visualization of hexagonal (left) and trigonal (right) 
perturbations to an ellipsoidal Fermi surface used to study 
the qualitative dependencies of the values of the rigidity co- 
efficients on the form of the Fermi surface. 

For k 3 > the upper critical field is described in the 
first order of perturbation theory by the second term in 
(p8|), and indeed the increase in K3 is accompanied by the 
appearance and growing of a "bump" around 9 = ir/2. 
The value H c2 (6 = §) equals (@). 

Non-zero coefficient K5 (trigonality) results in the dis- 
crimination between upper and lower parts of a crystal, 
and thus in the curves H C 2(6) non-symmetric with re- 
spect to the point 6 = 7r/2 if also the coefficient K3 7^ 
(see Fig. |). 

At non-zero K5 and K3 the curves are characterized 
by the presence of two maxima in the range < 8 < 
7r with the heights proportional to k 5 and the ratio of 
the heights proportional to |«3|, though at large |«a| the 
proportionality of the ratio becomes weaker, so that ( |2g| ) 
extrapolated to large |«a j exaggerates the difference in 
heights. 

A local minimum of H C 2{9) is achieved between the 
two maxima at an angle 8 < tt/2 for K3 < and 8 > ir/2 
for K3 > 0. 



V. RELATIVE VALUE OF THE RIGIDITY 
COEFFICIENTS Ki,k 2 ,k 3 ,K4 AND k 5 

In the weak=coupling BCS theory these coefficients can 
be calculated! 2 ] as averages over the Fermi surface 

» = (v Ft v Fj ^X k )^Xk))F, (32) 

Kq 

where V'u(^) are normalized vector basis functions (in the 
k-space) of the irreducible representation T according to 
which the order parameter 

dimT 

d i = E ^vvW (33) 

transforms. 

The details of the calculations are placed into Ap- 
pendix B. Here we describe the conclusions. 
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Table I: Rigidity coefficients for various basis functions and a spherical Fermi-surface with hexagonal and trigonal perturbations. 



Basis ipn(k) 



«2 _ "3 
K K 



(k x , k y 



III 

5 5 5 



U (U 2 -k 2 -9k k \v - 4- 256,^ _ 256 1 68 

l\iz\>^x ^y: ^f^xi^y 3 429 h cx 429 a ncx 3 429 tr 

£• fP - U Pi f£ A- 1z -^-4- -52-m, -5- 4- -5£-rn, - 2352„,, 

^2^3; Oh.x'^yJ \">x, i^y 26 4199 h cx 26 4199 hex 13 4199 hcx 221 tr 

i- (p - U Pi £ rP - P -2k A- 1z -2- - 4222 m. J. 1 4800 _5_ _ 4480 _ _42_ 

K> z yh, x Ohjxt^yJ ^z\^x ^yi ^f^xi^y 34 7429 hex 34 ' 7429 h cx 17 7429 h cx 323 tr 



The symmetry of a crystal can be taken into account 
in two ways. First, it is necessarily reflected in the sym- 
metry of the Fermi surface of a superconductor. And 
secondly, it should be complied with by the basis func- 
tions tp^ik) of the superconducting state. The shape of 
the Fermi surface (namely, its cross sections) is studied 
in the dHvA experiments. The basis functions are totally 
unobservable. 

We have modeled the Fermi surface by a sphere (for 
simplicity of calculation) with superimposed small hexag- 
onal and trigonal perturbations (see Fig. ^): 

/(*) = a he Jt(k x + iky) 6 + a t XMk + iky) 3 (34) 
= sin 3 9 (ahex sin 3 9 cos 6<p + an cos 9 sin 3y) , 

where ahex and atr are the small phenomenological pa- 
rameters defining the strength of the perturbation. 

The factor k z = cos 9 in the expression for the trigonal 
part of the perturbation differentiates between the up- 
per and lower halves of the Fermi surface. The form of 
the upper part of the perturbed surface is schematically 
shown with the solid curve on Fig. ^|, and that of the 
lower part — with the dashed one. The power of sin 9 in 
front of each term is of course a question of convenience, 
it must only be positive to make the two poles 9 = 0, it 
non-singular. 

We have taken four sample bases ipn(k) of the irre- 
ducible representation E u of the group D 3c t- The results 
are presented in Table |. The complicated fractional val- 
ues have no physical meaning, — e.g. for an oblate un- 
perturbed Fermi surface they would be different. Does 
matter the sensitivity of the rigidity coefficients of the su- 
perconducting state with a particular basis to the Fermi 
surface perturbation of a given symmetry. 

The first two bases are the most frequent examples for 
the E\ u and E^u ip-educible representations respectively 
of the group -D 6 JI3. Since both these representations 
become E M in the less symmetric group D^, their bases 
are valid choices. 

The second two bases are simply products of the first 
two by the function k z (k x — 3k x ky) = k z $l(k x + iky) 3 . 
Since this function is invariant under transformations 
only from the D^d group and not from D^h, these bases 



are specific for the E u representation of the D$d group. 

Table | shows that all four bases respond to trigonal 
perturbation of the Fermi surface, while only the two 
last respond to hexagonal perturbation, — basis (k x ,k y )z 
does not respond to the latter at all, and for basis k z (k\ — 
ky, —2k x k y )z the coefficient K4 rests unchanged in the 
linear order. 

VI. DISCUSSION 

As we have seen, six-fold modulations oc (1 — T/T c ) of 
the upper critical field in the basal plane can be naturally 
explained in the framework of a usual GL theory in a 
trigonal unconventional superconductor. Experimentally 
observed modulations (\SH^ X /H c2 \ < 0.02) are of the 
same order of magnitude as trigonal deviations from the 
hexagonal structure (~ 0.01). 

The form of the orientational dependence of H c2 as 
given by the theory rests fairly accurately sinusoidal in 
almost all the domain of the phenomenological parame- 
ters of the GL theory that keep the GL functional positive 
definite except for the values of both the parameters K3 
and K5 extremely close to the upper limit of their do- 
mains. 

At the same time, the curves of the angular depen- 
dence of resistivity having been observed in the transport 
experiments^ and their continuatiorJlJ, possessed sharp 
minima at <p = 7rn/3 (i.e., at directions of H || a). Not 
discussing here the legitimacy of mapping form of p(tp) 
onto H C 2((p) we remark here that the observed sharp dips 
can't be accounted for by the present theory of equilib- 
rium upper critical field. 

In contrast, the six- fold modulations in the theory 
developed by SaulsB in the framework of hexagonal 
structure of UPt3 originated from hexagonal in-plane 
anisotropy of the AFM order parameter m s coupled to 
the superconducting order parameter, on the supposition 
that m s stays almost orthogonal to the external magnetic 
field. 

Although the latter assumption contradicts experi- 
ment, the theory managed to account for the change of 
sign of the modulations with the passing through the 



7 




Figure 5: Schematic representation of the orientation of the 
coordinate system with respect to the crystalline axes. Also 
shown is the form (in the k-space) of the trigonal perturbation 
to the Fermi-surface. 



tetracritical point. Corrections from not exact orthogo- 
nality of m s to H also rendered hexagonal modulations 
non-sinusoidal. 

We considered a version of our GL theory with an 
orthorhombic anisotropy in Appendix A. Upper critical 
field retained six-fold modulations if the orthorhombic 
anisotropy was fixed in the direction parallel or perpen- 
dicular to the magnetic field, the sign of the modulations 
changing at the tetracritical point. However, the mod- 
ulations then became oc (1 — T/T c ) 2 , and the trigonal 
invariant oc K5 in the GL functional ([[]) then prevented 
the intersection of the two curves H C 2 (T) . 

Thus experimentally irrelevant hypothesis of m s stay- 
ing perpendicular to H leads in both theories to similar 
effects, although six- fold modulations of H C 2 themselves 
stem from completely different sources. 

Experimentstil revealed also correlation between the 
six-fold modulations of the upper critical field in the basal 
plane of some samples of UPt3 and the appearance of a 
non-symmetric peak at 9 = ir/2 (i.e., at H _L c) of H c2 
in the plane passing through the main symmetry axis. 

Since UPt3 turned out in reality to have trigonal crys- 
tal structure, one should expect both six-fold modula- 
tions of the upper critical field in the basal plane and 
some non-symmetry 9 «-> n — 9 in the plane of the sym- 
metry axis. But again, the form of the curves on Fig. || 
bears little similarity with usual hexagonal curves with 
superimposed non-symmetric peaks at 9 = ir/2 experi- 
mentally observed in the resistivity measurements o£J. 

An important and not yet fulUf understood role, how- 
ever, is played by the discovered!!] merohedry of UPt3 — 
division of a sample on a large number of small (~200- 
30fjA) domains with alternating sign of trigonality. 

The existence of such domains probably leaves the con- 
clusion of the modulation of H C 2 in the basal plane intact, 
for so six-fold a modulation is shared by all the domains 
irrespective of the signs of their trigonality. On the other 
hand, the magnitude of H c2 when a field is deflected from 
the basal plane clearly depends on the direction in which 



it is deflected — in other words, a trigonal monocrystal 
does not possess the symmetry c <-> — c. 

In the investigated samples total numbers of the do- 
mains with the opposite signs of trigonality roughly co- 
incided, so one should not probably wonder at the fact 
that the forms of the experimentally observedtil and the- 
oretically predicted (Fig. ^) dependencies are in poor 
agreement. Even the strong non-symmetry 9 <-> ir — 9 of 
the theoretical curves finds only modest correspondence 
in the experimental results probably also due to approx- 
imately equal number of merohedral domains with the 
opposite signs of trigonality. 

In the end, a few words are to be said on the micro- 
scopically obtained qualitative dependence of the trigonal 
GL coefficient Kg on the shape of the Fermi surface. We 
found that K5, which is responsible for all the effects of 
the trigonal structure in the theory, is proportional to a 
trigonal contribution to the model Fermi surface. The 
shape resembling right figure on Fig. f| was never ob- 
served to the author's knowledge in dHvA experiments 
in UPt3, probably because, as one can note, the cross- 
section of such a figure in linear order of the trigonal 
perturbation coincides with that of an unperturbed one. 
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Appendix A: VARIATIONAL METHOD FOR H C i 
IN THE BASAL PLANE 

Main features of the upper critical filed (|23| ) can be 
seen on a variational fit with the trial function 

= ( qiie- Vgg 1 / 2 \ 

71 I v±e -^Kl23/K4Z 2 /2 V ' > 

where r}± and r]\\ are to be found from minimization of 
the GL functional with respect to them. The result then 
turns out to be just a truncated sum of the second order 
perturbation theory. 

Even simpler though for this particular form of the trial 
function istl to substitute it directly into Eq. (0) , which 
then becomes an algebraic equation on rj±, m. From the 
condition of vanishing of the corresponding determinant 
one then gets for the eigenvalue 

(A - ^KiK 4 )(A - V K 123«4) = V K 1 K 4V K 123K4^4 2 , (A2) 
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where 

A 2 = 



2k 2 



^'^^ cos^. (A3) 



K4 ^/KiKi2 3 (y/Kl + V«123) 

Using the definition ([?]) of A, one arrives at an equation 

(H - ff c2 ||(T)) (H - H c2X (T)) = A 2 H 2 , (A4) 

which describes six-fold modulation of the upper critical 
field. E.g., supposed that K23 > 0, i.e., that the super- 
conducting state with the order parameter parallel to the 
magnetic field arises first, we find 



H, 



c2 



H, 



c2\\ 



1 



A 2 H, 



c2\\ 



H, 



c2|| 



H, 



c2± 



(A5) 



or, substituting (jlj) and (fL5|), 

H C 2 



H, 



''211 



1 + 2^ 



«123 V K 123 



Kl (V«123 + V^l) 



COS 3cO 



(A6) 



For K23 < we get an analogous expression with the 
opposite sign in front of the corrections. 

It is interesting to remark that a six-fold modulation 
of H C 2 in the basal plane of a trigonal crystal appears 
alsoja|-the presence of an orthorhombic anisotropy (see, 
e.g.,tll,E3) in the direction fixed parallel or perpendicular 
to the magnetic field. The corresponding term was in- 
troduced to the GL theory to account for the splitting of 
the superconducting transition in UPt3 . 

With such a term one of the fields H c2 \\ (T) or H c2 ±(T) 
becomes oc (1 — T/T c «) rather than oc (1 — T/T c ), where 
the splitting between T c and T„ is proportional to the 
magnitude of the orthorhombic anisotropy. Eq. (A4) 
is still valid then and produces six-fold corrections (AE) 
to H c2 , though proportional to (1 — T/T c ) 2 unlike (Af), 
where they are linear. In such a version of the theory, 
however, the term oc K5 in the GL functional leads to the 
repulsion of the two lines of H C 2(T) so that instead of a 
tetracritical point one has a splitting oc A. 



Appendix B: CALCULATION OF THE RIGIDITY 
COEFFICIENTS 

Relative values of the rigidity coefficients is convenient 
to study with the model dispersion law of electrons 



e k = ~z h 2e F /(fc), 
Am 



(Bl) 



where f(k) is a small perturbation with a symmetry de- 
termined by the symmetry of the crystal. For simplicity 
we will consider all the expressions up to leading term in 
/. The form of the perturbation of interest to us is (|3^ ) 
as stated in the main text. 

The Fermi-surface is defined by £k = £f, whence we 
obtain its equation in the parameterization of the unit 
vector k 



k(*o = Mi - /(*))*, 



(B2) 



where &f = ^/2mep. Making use of the value of k in 
terms of the spherical angles 9, if one obtains for the el- 
ement of the Fermi surface 



|<9 e k x d v k\d9dtp = jfc|(l - 2 f(k)) sin BdBdip. (B3) 



The unit vector vp of the Fermi velocity vp = SkCklp 
can be calculated from ([Bl]) to be 



Vf 



k + d~j-k(kd % j) 



(B4) 



Thus the averaging ( |32[) over the Fermi-surface may be 
written as 



((1 - 2f{k))v Fl v Fo i>^k)xl, v {k)) k 



(B5) 



(l-2/(fe)) S 
= (l + 2</(%) {kik^^Mk 

+ (ft(V-*i(*V)-*i/) 

+ kihj - Hkd k f) - fcj)]^(fc)^(fc) ) 



Calculating the integrals in (B5) with the bases func- 
tions described in the main text, one arrives at Table §. 
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